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QCD Description of Hadronic Tau Decays 

A. Pich^ 

'"Departament de Fi'sica Teorica, IFIC, Univ. Valencia-CSIC, Apt. 22085, E-46071 Valencia, Spain 

The QCD analysis of hadronic r decays is reviewed and a summary of the present phenomenological status is 
presented. The following topics are discussed: the determination of as{m'^) = 0.338 ± 0.012 from the inclusive r 
hadronic width, the measurement of |Kis| through the Cabibbo-suppressed decays of the r, and the extraction of 
chiral-perturbation-theory couplings from the spectral tau data. 



1. THEORETICAL FRAMEWORK 

The hadronic r decays turn out to be a beau- 
tiful laboratory for studying strong interaction 
effects at low energies [1-3 . The r is the 
only known lepton massive enough to decay into 
hadrons. Its semileptonic decays are then ideally 
suited to investigate the hadronic weak currents. 

The inclusive character of the total r hadronic 
width renders possible an accurate calculation of 
the ratio gHS] 



Rr — 



r [t Vt- hadrons] 



r[T- 



Rt,V + Rt.A + Rt.S 



The theoretical analysis involves the two-point 
correlation functions for the vector — '4!jj^ipi 
and axial-vector — '4'jJ^'y5ipi colour-singlet 
quark currents — u, d, s): 



which have the Lorentz decompositions 



(2) 



where the superscript (J = 0,1) denotes the an- 
gular momentum in the hadronic rest frame. 



The 



imagmary parts of n[^|y(g^) are pro- 
portional to the spectral functions for hadrons 
with the corresponding quantum numbers. The 
hadronic decay rate of the r can be written as 
an integral of these spectral functions over the 



invariant mass s of the final-state hadrons: 



ds 



Rr ^ 12tt 1 - 
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The appropriate combinations of correlators are 
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The contributions coming from the first two terms 
correspond to Rry i?r,A respectively, while 
Rt,s contains the remaining Cabibbo-suppressed 
contributions. 

The integrand in Eq. ([3]) cannot be calculated 
at present from QCD. Nevertheless the integral it- 
self can be calculated systematically by exploiting 
the analytic properties of the correlators n'^''^(s). 
They are analytic functions of s except along the 
positive real s-axis, where their imaginary parts 
have discontinuities. Rt can then be written as a 
contour integral in the complex s-plane running 

= ml [5]: 



counter-clockwise around the circle 
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correlators only 
which is signifi- 



This expression requires 
for complex s of order 
cantly larger than the scale associated with non- 
perturbative effects. Using the Operator Product 
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Expansion (OPE), n('')(s) = Ez? C*!/ 
to evaluate the contour integral, Rt can be ex- 
pressed as an expansion in powers of The 
uncertainties associated with the use of the OPE 
near the time-like axis are heavily suppressed by 
the presence in ([5|) of a double zero at s = m^. 

In the chiral limit {mu,d,s = 0), the vector 
and axial-vector currents are conserved. This 
implies sn'^''^(s) = 0. Therefore, only the cor- 
relator n(°+i)(s) contributes to Eq. bmce 
(1 - xfil + 2x) = 1 - 3x2 -f 22-3 = s/ml], 
Cauchy's theorem guarantees that, up to tiny 
logarithmic running corrections, the only non- 
perturbative contributions to the circle integra- 
tion in ([5]) originate from operators of dimensions 
D = 6 and 8. The usually leading D = 4 opera- 
tors can only contribute to Rt with an additional 
suppression factor of 0(q;^), which makes their 
effect negligible [B]. 

2. DETERMINATION OF 

The Cabibbo-allowed combination Rry+A can 
be written as [6] 

Rr.V+A = Nc iKdP ^EW {1 + (Jp + <5np} , (6) 

where Nc = 3 is the number of quark colours and 
Sew = 1-0201 ± 0.0003 contains the electroweak 
radiative corrections [9HTT|. The dominant cor- 
rection 20%) is the perturbative QCD con- 
tribution 6p, which is already known to O(ck^) 
[HI [12]. Quark mass effects [B1[T31[I1] are tiny 
for the Cabibbo-allowed current and amount to 
a negligible correction smaller than 10"'' [HjllSj. 

Non-perturbative contributions are suppressed 
by six powers of the r mass [6] and, therefore, 
are very small. Their numerical size has been 
determined from the invariant-mass distribution 
of the final hadrons in r decay, through the study 
of weighted integrals [TBI . 

which can be calculated theoretically in the same 
way as Rr, but are more sensitive to OPE cor- 
rections. The predicted suppression (Bl of the 
non-perturbative corrections to Rr has been con- 



firmed by ALEPH [H], CLEO [H] and OPAL 
[T^. The most recent analysis gives [3D] 

(5np = -0.0059 ±0.0014. (8) 

The QCD prediction for Rt,v+a is then com- 
pletely dominated by Sp ; non-perturbative effects 
being smaller than the perturbative uncertain- 
ties from uncalculated higher-order corrections. 
Assuming lepton universality, the measured val- 
ues of the r lifetime and leptonic branching ra- 
tios imply Rr = 3.6291 ± 0.0086 [H]. Subtract- 
ing the Cabibbo-suppressed contribution i?r,s = 
0.1613 ± 0.0028 21j, one obtains Rr.v+A = 
3.4678 ±0.0090. Using |K<i| = 0.97425 ± 0.00022 
[22j and ([8]), the pure perturbative contribution 
to Rr is determined to be: 

(5p = 0.1997 ±0.0035. (9) 

The predicted value of 5p turns out to be very 
sensitive to Q;s(m^), allowing for an accurate de- 
termination of the fundamental QCD coupling 
[51[5j. The calculation of the ©(a*) contribu- 
tion [12j has triggered a renewed theoretical in- 
terest on the asim^) determination, since it al- 
lows to push the accuracy to the four-loop level. 
However, as shown in Table [1] the recent theoret- 
ical analyses slightly disagree on the final result. 
The differences are larger than the claimed 0{a^) 
accuracy and originate in the different inputs or 
theoretical procedures which have been adopted. 

2.1. Perturbative contribution to Rr 

In the chiral limit, the result is more con- 
veniently expressed in terms of the logarithmic 
derivative of the two-point correlation function 
of the vector (axial) current, n(s) = i n*^°+'-'(s), 
which satisfies an homogeneous renormalization- 
group equation: 

n=0 ^ ^ 

With the choice of renormalization scale ^ —s 
all logarithmic corrections, proportional to pow- 
ers of log(— s/z^t^), have been summed into the 
running coupling. For three flavours, the known 
coefficients take the values: Kg = Ki = 1; K2 = 
1.63982; KaCMS) = 6.37101 and KiiMS) = 
49.07570 [ig. 
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Table 1 



C'(af) determinations of as{m'^). The assumed value of dp is also given (if quoted by the authors). 



Reference 


Method 


Sp 






Baikov et al. [H] 


CIPT, FOPT 


0.1998 ±0.0043 


0.332 ±0.016 


0.1202 ±0.0019 


Davier et al. [53] 


CIPT 


0.2066 ±0.0070 


0.344 ± 0.009 


0.1212 ±0.0011 


Beneke-Jamin [13] 


BSR + FOPT 


0.2042 ± 0.0050 


0.316 ±0.006 


0.1180 ± 0.0008 


Maltman-Yavin [24] 


PWM + CIPT 




0.321 ±0.013 


0.1187 ±0.0016 


Menke [25] 


CIPT, FOPT 


0.2042 ± 0.0050 




0.1213 ±0.0012 


Caprini-Fischer [26] 


BSR + CIPT 


0.2042 ± 0.0050 


320 + ° °" 




Cvetic et al. [27] 


/3exp + CIPT 


0.2040 ±0.0040 


0.341 ± 0.008 


0.1211 ±0.0010 


Pich [Tj 


CIPT 


0.2038 ±0.0040 


0.342 ±0.012 


0.1213 ±0.0014 



Table 2 

Exact results for A*^")(as) {n < 4) at different /3- function approximations, and corresponding values of 
6p = J2t=i A^'^Kots), for ar = as{ml)/TT = 0.11. The last row shows the FOPT estimates at O(a^). 







v4W(a,) 








Sp 


Pn>l — 





0.14828 


0.01925 


0.00225 


0.00024 


0.20578 


Pn>2 — 





0.15103 


0.01905 


0.00209 


0.00020 


0.20537 


/3n>3 — 





0.15093 


0.01882 


0.00202 


0.00019 


0.20389 


/3n>4 — 





0.15058 


0.01865 


0.00198 


0.00018 


0.20273 


o{4) 




0.16115 


0.02431 


0.00290 


0.00015 


0.22665 



The perturbative component of Rr is given by 
<5p - ^if„A(")(a.), (11) 

n=l 

where the functions [7] 
A(")(a.).-L/ dsra^Y 

2tTI J\s\=rnl S \ TT J 

X fl-24 + 24-4) (12) 

are contour integrals in the complex plane, which 
only depend on Cr = as{m'^)/Tr. Using the 
exact solution (up to unknown /3„>4 contribu- 
tions) for as{—s) given by the renormalization- 
group /3-function equation, they can be numer- 
ically computed with a very high accuracy [3- 
Table [2] gives the numerical values for A^"^as) 
(n < 4) obtained at the one-, two-, three- and 
four-loop approximations (i.e. /3„>i = 0, /3„>2 = 
0, l3n>3 — and /3n>4 ~ 0, respectively), to- 
gether with the corresponding results for Sp = 



Y,t=i if„^("Has), taking a, = 0.11. The per- 
turbative convergence is very good and the results 
are stable under changes of the renormalization 
scale. The error induced by the truncation of the 
[3 function at fourth order can be conservatively 
estimated through the variation of the results at 
five loops, assuming = ±(31/(33 — t443, i.e. a 
geometric growth of the f3 function. 

Higher-order contributions to the Adler func- 
tion D{s) will be taken into account adding the 
fifth-order term K5 A^^) (a^) with K5 = 275±400. 
Moreover, we will include the 5-loop variation 
with changes of the renormalization scale in the 
range /z^/(— s) G [0.5,1.5]. Adopting this very 
conservative procedure, the experimental value of 
Sp given in Eq. ^ implies 

as(TO^) = 0.338 ±0.012. (13) 

The result is slightly lower than the one given in 
Ref. [1], due to the smaller value of Sp. 

The strong coupling measured at the t mass 
scale is significantly larger than the values ob- 
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tained at higher energies. From the hadronic de- 
cays of the Z, one gets a,(M|) = 0.1190 ±0.0027 
[25] . which differs from as{m'^) by 18a. After 
evolution up to the scale Mz [29] , the strong cou- 
pling constant in (IT3t decreases to 

as(Af|) = 0.1209 ±0.0014, (14) 

in excellent agreement with the direct measure- 
ments at the Z peak and with a better accuracy. 
The comparison of these two determinations of as 
in two very different energy regimes, rrir and Mz, 
provides a beautiful test of the predicted running 
of the QCD coupling; i.e., a very significant ex- 
perimental verification of asymptotic freedom. 

2.2. Fixed-order perturbation theory 

The integrals A'-'^^as) can be expanded in 
powers of a^, ^("H^s) < + C(a?+^)- One 
recovers in this way the naive perturbative ex- 
pansion [7] 

n— 1 n—1 

This approximation is known as fixed-order per- 
turbation theory (FOPT), while the improved ex- 
pression (fTTj) . keeping the non-expanded values of 
v4*^"^(q!s), is usually called contour-improved per- 
turbation theory (CIPT) [7l|3^. 

As shown in the last row of Table [21 even at 
0{a'^), FOPT gives a rather bad approximation 
to the integrals A*^")(q:s), overestimating Sp by 
12% at Or = 0.11. The long running of as{—s) 
along the circle \s\ = ra^ generates very large gn 
coefficients, which depend on Km<n and Pm<n 
[7]: gi = 0, 52 = 3.56, 53 = 19.99, 54 = 78.00, 
35 = 307.78. These corrections are much larger 
than the original Kn contributions and lead to 
values of as{rn^) smaller than ([T3|l . FOPT suffers 
from a large renormalization-scale dependence [7] 
and its actual uncertainties are much larger than 
usually estimated [25]. 

The origin of this bad behaviour can be un- 
derstood analytically at one loop [7]. In FOPT 
one makes within the contour integral the series 
expansion (log (—s/m^) = i(j), (j) ^ [— tt, tt]) 



which is only convergent for Ot- < 0.14. At 
the four-loop level the radius of convergence is 
slightly smaller than the physical value of a^- 
Thus, FOPT gives rise to a pathological non- 
convergent series. The long running along the 
circle makes compulsory to resum the large loga- 
rithms, log" (— s/rn^), using the renormalization 
group. This is precisely what CIPT does. 

2.3. Renormalon hypothesis 

The perturbative expansion of the Adler func- 
tion is expected to be an asymptotic series. If its 
Borel transform, B{t) = X]n=o ^"-1-1*"/"'' '^ere 
well-behaved, one could define D{s) through the 
Borel integral 

m = |l + dt e--*/"=(^)i?(i)} . (17) 

However, B{t) has pole singularities at positive 
(infrared renormalons) and negative (ultravio- 
let renormalons) integer values of the variable 
u = — /3it/2, with the exception of it = 1 |31) . 
The infrared renormalons at m = +n are related 
to OPE corrections of dimension D = 2n. The 
renormalon poles closer to the origin dominate 
the large-order behaviour of D{s). 

It has been argued that, once in the asymptotic 
regime (large n), the rcnormalonic behaviour 
of the Kn coefhcients could induce cancelations 
with the running gn corrections, which would be 
missed by CIPT. In that case, FOPT could ap- 
proach faster the 'true' result provided by the 
Borel summation of the full renormalon series 
(BSR) [15]. This happens actually in the large- 
/?! limit :32, 33 , which however does not approxi- 
mate well the known Kn coefficients. A model of 
higher-order corrections with this behaviour has 
been recently advocated [15]. The model mixes 
three different types of renormalons {n = —1, 2 
and 3) plus a linear polynomial. It contains 5 free 
parameters which are determined by the known 
values of -fs'1.2,3,4 and the assumption K^ = 283. 
One gets in this way a larger ^p, implying a 
smaller value for as(m^). The result looks how- 
ever model dependent [34] . 

The implications of a rcnormalonic behaviour 
have been put on more solid grounds, using an 
optimal conformal mapping in the Borel plane. 
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which achieves the best asymptotic rate of conver- 
gence, and properly implementing the CIPT pro- 
cedure within the Borel transform |26| . Assuming 
that the known fourth-order series is already gov- 
erned by the u = —1 and u = 2 renormalons, the 
conformal mapping generates a full series expan- 
sion (i^5 = 256, Ke = 2929 . . . ) which results, 
after Borel summation, in a larger value of Sp; i.e. 
the Kn>4 terms give a positive contribution to Sp 
implying a smaller asim^) [26j . 

Renormalons provide an interesting guide to 
possible higher-order corrections, making appar- 
ent that the associated uncertainties have to be 
carefully estimated. However, one should keep in 
mind the adopted assumptions. In fact, there are 
no visible signs of renormalonic behaviour in the 
presently known series: the n — —1 ultraviolet 
renormalon is expected to dominate the asymp- 
totic regime, implying an alternating series, while 
all known Kn coefhcients have the same sign. One 
could either assume that renormalons only be- 
come relevant at higher orders, for instance at 
n = 7, and apply the conformal mapping with 
arbitrary input values for and Kq. Differ- 
ent assumptions about these two unknown coefh- 
cients would result in different central values for 
as{ml). 

A different reshuffling of the perturbative se- 
ries, not related to renormalons, has been recently 
proposed [27] . Instead of the usual expansion in 
powers of the strong coupling, one expands in 
terms of the /3 function and its derivatives (/3cxp), 
which effectively results in a different estimate of 
higher-order corrections. One gets in this way a 
weaker dependence on the renormalization scale 
and a value of as{m'^) similar to the standard 
CIPT result. 

2.4. Non-perturbative corrections 

At the presently achieved precision, one should 
worry about the small non-perturbative correc- 
tions. In fact, a proper definition of the infrared 
renormalon contributions is linked to the corre- 
sponding OPE corrections with D — 2n. A re- 
cent re-analysis of the ALEPH data [24,, with 
pinched-weight moments of the hadronic distri- 
bution (PWM) and CIPT, obtains ^^(m^) = 
0.321 ± 0.013. This smaller value originates in 



a different estimate of the non-perturbative con- 
tributions. Unfortunately, Ref. [Mj does not 
quote any explicit values for (5np and Sp. From 
the information given in that reference, I deduce 
Jnp = 0.012 ± 0.018. Ahhough compatible with 
([8]), the central value is larger and has the op- 
posite sign. This shift implies a smaller Sp and, 
therefore, a slightly smaller strong coupling. 

The so-called duality violation effects, i.e. the 
uncertainties associated with the use of the OPE 
to approximate the exact correlator, have been 
also investigated [23l|35j. Owing to the presence 



in ([5]) of a double zero at s = m^, these effects 
are quite suppressed in Rr . They are smaller than 
the errors induced by S^p , which are in turn sub- 
dominant with respect to the leading perturbative 
uncertainties. 



3. |Vus DETERMINATION 

The separate measurement of the \AS\ = 
and I AS* I = 1 tau decay widths provides a very 
clean determination of Vus [361137] . To a first ap- 
proximation the Cabibbo mixing can be directly 
obtained from experimental measurements, with- 
out any theoretical input. Neglecting the small 
SU(3)-breaking corrections from the — rrid 
quark-mass difference, one gets: 



Rt.. 



Rt,V+A 



1/2 



0.210±0.002. 



The new branching ratios measured by BaBar 
and Belle are all smaller than the previous world 
averages, which translates into a smaller value 
of Rt,s and |Kis|- For comparison, the previ- 
ous value Rr,s = 0.1686 ± 0.0047 20 resulted 
in IKsP"^^^ = 0.215 ±0.003. 

This rather remarkable determination is only 
slightly shifted by the small SU(3)-breaking con- 
tributions induced by the strange quark mass. 
These effects can be estimated through a QCD 
analysis of the differences [ni[TI[ l55W5] 



SRt' 



r>kl 



Rr,S 
I ''lis 



(18) 



The only non-zero contributions are proportional 
to the mass-squared difference rrij. — or to 
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vacuum expectation values of SU(3)-breaking op- 
erators such as SO4 = {0\msSS — mddd\Q) k, 
(-1.4 ± 0.4) • 10-3 GeV" The dimen- 

sions of these operators are compensated by cor- 
responding powers of m^, which imphes a strong 
suppression of 5R^^ [13] : 

5Rf « 245ew (l^iM) 

-27r2^g,,(a.)| , (19) 

where = md/m^ = 0.053 ± 0.002 iM]- The 
perturbative corrections Afei(as) and Qfei(as) are 
known to 0(0^) and ©(a^), respectively p ^ fM ] . 

The J = contribution to Aoo(q:s) shows a 
rather pathological behaviour, with clear signs of 
being a non-convergent perturbative series. For- 
tunately, the corresponding longitudinal contri- 
bution to 5Rr = 5R^^ can be estimated phe- 
nomenologically with a much better accuracy, 
5Rr\^ = 0.1544 ± 0.0037 [3i|45], because it is 
dominated by far by the well-known r — >■ Vt-t: 
and r — VtK contributions. To estimate the re- 
maining transverse component, one needs an in- 
put value for the strange quark mass. Taking the 
range m,{mr) = (100 ± 10) MeV [ms(2 GeV) = 
(96 ± 10) MeV], which includes the most recent 
determinations of from QCD sum rules and 
lattice QCD gg, one gets finally SRr,th = 0.216± 
0.016 |37], which implies 

1^ I - ^ ^ 

= 0.2166 ±0.0019oxp± 0.0005th. (20) 

A larger central value, |Ks| = 0.2217 ± 0.0032, is 
obtained with the old world average for i?r,s- 

Sizeable changes on the experimental determi- 
nation of i?T,s could be expected from the full 
analysis of the huge BaBar and Belle data sam- 
ples. In particular, the high-multiplicity decay 
modes are not well known at present. The re- 
cent decrease of several experimental tau branch- 
ing ratios is also worrisome. As pointed out by 
the PDG [22], 15 of the 16 branching fractions 
measured at the B factories are smaller than the 
previous non-B-factory values. The average nor- 



malized difference between the two sets of mea- 
surements is —1.36(7. Thus, the result ([20]) could 
easily fluctuate in the near future. In fact, com- 
bining the measured Cabibbo-suppressed r dis- 
tribution with electroproduction data, a slightly 
larger value of \Vus\ is obtained [46] . 

The final error of the Vus determination from 
T decay is dominated by the experimental uncer- 
tainties. If Rt.s is measured with a 1% precision, 
the resulting V^g uncertainty will get reduced to 
around 0.6%, i.e. ±0.0013, making r decay the 
best source of information about Vus- 

An accurate measurement of the invariant- 
mass distribution of the final hadrons could make 
possible a simultaneous determination of Vus 
and the strange quark mass, through a corre- 
lated analysis of several weighted differences SRf . 
However, the extraction of suffers from theo- 
retical uncertainties related to the convergence of 
the perturbative series Aki{as). A better under- 
standing of these corrections is needed. 

4. CHIRAL SUM RULES 

When mu^d.s — 0, the QCD Lagrangian has 
an independent SU{3) flavour invariance for the 
left and right quark quiralities. The two quirali- 
ties have exactly the same strong interaction, but 
they are completely decoupled. This chiral in- 
variance guarantees that the two-point correla- 
tion function of a left-handed and a right-handed 
quark currents, Ulr{s) = n|f+^^(s) - n^°d+j'(s), 
vanishes identically to all orders in perturbation 
theory (the vector and axial- vector correlators re- 
ceive identical perturbative contributions). The 
non-zero value of IlLRis) originates in the spon- 
taneous breaking of chiral symmetry by the QCD 
vacuum. At large momenta, the corresponding 
OPE only receives contributions from operators 
with dimension d > 6, 




The non-zero up and down quark masses induce 
tiny corrections with dimensions two and four, 
which are negligible at high energies. 

At very low momenta, Chiral Perturbation 
Theory (xPT) dictates the low-energy expansion 
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ofIlLj^{s) in terms of the pion decay constant and 
the xPT coupUngs [Oip*)] and Cs7 

Analyticity relates the short- and long-distance 
regimes through the dispersion relation 

-— / dsw{s)ULR{s) = - / dsw{s)p{s) 
+ 2/2 w{ml) + Res[w{s)nLBis),s = 0] , (22) 

where p{s) = ^ImllLRis) and w{s) is an arbi- 
trary weight function that is analytic in the whole 
complex plane except in the origin (where it can 
have poles). The last term in (|22p accounts for 
the possible residue at the origin. 

For So < m^, the integral along the real axis 
can be evaluated with the measured tau spectral 
functions. Taking w{s) = s" with n > 0, there is 
no residue at the origin and, with sq large enough 
so that the OPE can be applied in the entire circle 
\s\ = So, the OPE coefficients are directly related 
to the spectral function integration. With n = 
and 1, there is no OPE contribution in the chiral 
limit and one gets the celebrated first and second 
Weinberg sum rules [17]. For negative values of 
the integer n, the OPE does not contribute ei- 
ther while the residues at zero are determined by 
the xPT low-energy couplings, which can be then 
experimentally determined [48] . 

Moreover, the absence of perturbative contri- 
butions makes ((22)) and ideal tool to investigate 
possible quark-hadron duality effects, formally 
defined through PHST] 

BV^ = ^ (f ds w{s) (ULRis) - ngPE(5)) 

POO 

= / ds w{s) p{s) . (23) 

J Bo 

This has been thoroughly studied in Ref . [S2] , us- 
ing for the spectral function beyond Sz ~ 2.1 GeV 
the parametrization [551H51 - I5T| 

p(s>s,) = Ke-^^sin(/3(s-s,)), (24) 

and finding the region in the 4-dimensional 
{k,^, l3, Sz) parameter space that is compatible 
with the most recent experimental data |17| and 
the following theoretical constraints at so — > oo: 
first and second Weinberg sum sules [17] and the 



sum rule of Das et al. [53] that determines the 
pion electromagnetic mass difference. 

Ref. [52] performs a statistical analysis, scan- 
ning the parameter space (k,7, /3, Sz) and select- 
ing those 'acceptable' spectral functions which 
satisfy the experimental and theoretical con- 
straints. From a generated initial sample of 
160,000 tuples, one finds 1, 789 acceptable distri- 
butions compatible with QCD and the data. The 
differences among them determine how much free- 
dom is left for the behaviour of the spectral func- 
tion beyond the kinematical end of the t data. 
For each acceptable spectral function one calcu- 
lates the parameters Liq, Cgy, Og and Og, ob- 
tained through the dispersion relation with 
the appropriate weight functions. The resulting 
statistical distributions determine their finally es- 
timated values; the dispersion of the numerical 
results provides a good quantitative assessment 
of the actual uncertainties. 

The study has been also performed with 
pinched weight functions of the form w{s) = 
s"(s — Sz)™ (m > 0) that vanish at s = s^. 
As expected, these weights are found to mini- 
mize the uncertainties from duality-violation ef- 
fects, allowing for a more precise determination of 
the hadronic parameters. One finally obtains [52] , 

L^o(Afp) = -(4.06 ± 0.39) • 10"\ 

Q7(A/p) = (4.89 ± 0.19) ■ 10-^ GeY-^ , 

e>6 - (-4.3 1 J] ?) • 10^3 Qgv^ 

Os = (-7.2 ifj) • 10~3 GeV* . (25) 

The determination of the two xPT couplings is 
in good agreement with (but more precise than) 
recent theoretical calculations, using Resonance 
Chiral Effective Theory [S4^ and large-A''c tech- 
niques at the next-to- leading order [55], which 
predict [56]: ilo(Mp) = -(4.4 ± 0.9) • lO'^ and 
Clj{Mp) = (3.6±1.3)-10-3 GeV~^ It also agrees 
with the present lattice estimates oiL\Q{Mp) [57] . 

Duality- violation effects have very little impact 
on the determination of Lio and Cg? because the 
corresponding sum rules are dominated by the 
low-energy region where the data sits. Thus, one 
obtains basically the same results with pinched 
and non-pinched weight functions. This is no- 
longer true for Og and Og, which are sensitive to 
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Figure 1. Published results for Oe and Os [52] . 



the high-energy behaviour of the spectral func- 
tion; pinched-weights provide then a much better 
accuracy. This could explain the numerical dif- 
ferences among previous estimates pTlfTOlHSllSB] , 
shown in Fig. [1] where duality violation uncer- 
tainties were not properly assessed. The results 
fix with accuracy the value of Oq and deter- 
mine the sign of Os- This information is needed 
to calculate the electromagnetic penguin contri- 
bution to the CP-violating ratio e'^Iek [59] . 
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